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In this paper we investigate circular optical resonators with gaussian 
profiles of the mirror reflectivities. Closed form solution to the integral 
equations for such resonators are obtained. The dominant TEM 0i0 mode 
characteristics of a resonator consisting of one variable reflectivity mirror 
(VRM) and one unifoiin reflectivity mirror (URM) are considered in 
detail for a variety of parameters. This resonator is particularly suitable 
for high-gain lasers. Its advantages in comparison to the conventional type 
are: (i) there is larger mode volume utilization, and (ii) the power trans- 
mitted at the variable reflectivity mirror can in principle be utilized as the 
power output. We discuss dependence of the spot sizes on laser gain and 
mirror-curvature tolerances and present a specific design of a Fabry-Perot 
resonator for fundamental mode operation and the expected •performance. 

I. INTRODUCTION 

The dominance of the fundamental mode in optical resonators with 
uniform reflectivity mirrors is due to the lowest diffraction loss of this 
mode. The power output of this mode is commonly obtained by using 
a partially transparent mirror. These two features could be combined 
in a resonator consisting of one uniform reflectivity mirror (URM) 
and one variable reflectivity mirror (VRM). 

Resonators with VRM were investigated previously. S. N. Vlasov 
and V. I. Talanov 1 considered symmetrical two-dimensional resona- 
tors with two types of variations of the mirror reflectivities including 
the gaussian and obtained solutions for the eigenvalues of the resonator 
integral equations. N. G. Vakhimov 2 investigated the natural resonant 
frequencies and field distributions of symmetrical resonators with 
gaussian VRM by using an asymptotic method of solution to the wave 
equation subject to impedance boundary conditions. N. Kumagai and 
others 3 investigated Fabry Perot resonators with VRM of finite dimen- 
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sions by solving numerically the resonator integral equation for differ- 
ent mirror reflectivities. Y. Suematsu and others 4 studied beam wave- 
guides with gaussian transmission filters for the improvement of the 
stability of beam transmission. 

In this work we investigate nonsymmetrical circular resonators con- 
sisting of one URM and one VRM. The radii of curvature of the 
mirrors are arbitrary. The reflection coefficients of the VRM are as- 
sumed to have gaussian profiles in the radial direction. For such 
resonators with infinite mirrors, solutions to the resonator integral 
equations are obtained in terms of Laguerre functions with complex 
arguments. The modal fields decrease off-axis very rapidly and con- 
sequently these solutions are also applicable to resonators with finite 
mirrors. 

Resonators of the type considered seem to be particularly suitable 
for high-gain lasers as for example C0 2 lasers. It is shown subsequently 
that for the fundamental TE n .o mode, the spot sizes obtainable are 
considerably larger than those obtained with URM resonators of the 
same length and the same fundamental mode threshold gain ratio. 
This should result in a larger mode volume utilization. Furthermore, 
the power loss due to the transparency of VRM could be utilized as 
the power output. 

In the following sections the solutions for the resonator modes and 
eigenvalues of nonsymmetrical resonators are obtained. It is shown that 
the solutions for symmetrical resonators consisting of two identical 
VRM are readily obtainable as special cases. Mode-stability criteria 
are established as functions of the resonator geometries and VRM 
parameters. The spot sizes of the fundamental TE , mode are com- 
puted as a function of the threshold-gain ratio for a variety of pa- 
rameters. A comparison is made between the obtainable spot sizes 
with Fabry Perot resonators with VRM and URM. We show that 
much larger spot sizes can be achieved with the VRM resonator. 
We show also that the spot-size diameters of VRM mirrors depends 
basically on the threshold-gain ratio and on the mirror-curvature tol- 
erances. A specific design of a Fabry Perot resonator with a VRM is 
examined and the expected performance is presented. 

II. NONSYMMETRICAL RESONATORS 

2.1 Solutions to the Integral Equation (IE) 

The geometry of the nonsymmetrical resonator is shown in Fig. 1. 
It consists of one mirror with variable reflectivity (Ml) and one mirror 
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MIRROR 1 (Ml) MIRR0R2 (M2) 

r = r exp(-£ P 2 ) r 2 =i,o 

Fig. 1 — Nonsymmetrical resonator. 

with uniform reflectivity (M2). The separation between the mirrors 
is d, and the radii of curvature are designated by R x and R 2 respectively. 
The reflection coefficients of the VRM, T, is assumed to vary in the 
radial direction p as follows: 



r = r exp (-/3 P 2 ) 



(1) 



where T and /3 are constants with | r | ^ 1. 

The reflection coeflicient of the URM is assumed to be unity. (A 
reflection coeflicient different than unity can readily be included in the 
solution.) 

The integral equations for this resonator are obtained in a manner 
analogous to a URM resonator, by imposing the condition that the 
field should reproduce itself after a round trip. With the azimuthal 
dependence for the electric field E(p, <p) = exp (—jrf>)F t (p), the two 
simultaneous integral equations are: 5 

K W } (pO = i' +1 exp (-jkd)M f *T(pi) exp [-jM{g lP \ + g 2 p\)/2] 

■J t (M Pl p 2 )p 1 d Pl , (2) 

K W(pO = j M exp (-jkd)MT exp (-flp?) f F%p 2 ) 

Jo 

■ exp [-jM{g lP \ + ^ 2 p 2 2 )/2] • J t (M Pl p 2 )p 2 dp 2 , (3) 



2352 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1970 

where F^ipx), f} a, (p») are the radial field distributions at (Ml) and 
(M2) respectively, K™ and K ( t 2) are the associated eigenvalues, J t 
is a Bessel function of order I, M = 2ir/\d, X is the free-space wave- 
length and 

"-fr-D- (4) 

"-('-It)- (5) 

The integral equations (2) and (3) are solved by using the self- 
reciprocal properties of the Laguerre functions on the Hankel trans- 
form. 6,7 These properties are 

[ x v+1 exp (-0x 2 )L V n(ax 2 )J,(xy)(xy) h dx 

Jo 

= 2— V + *C8 - aW ' exp (-y 2 /^)^ ^^ J (6) 

where Z^ is a Laguerre function of order ^ w. 

Based on equation (6), the modal solutions to the integral equations 
are: 

F?\ Pl ) = exp (-y aP ?/2)I<»(«iP?)(V^ P.)', (7) 

F\ 2 \p 2 ) - exp (-7.pI/2)Lb(aipJ)(V^ /»)'■ (8) 

After some algebraic manipulations, the following relations are ob- 
tained for the parameters. 

Yi - «i + 0, (9) 

. M 2 . 



4?2 



^ + ^^-i/j 



(10) 



It is convenient to express «i in terms of a complex trigonometric func- 
tion as follows 

M 

a x = x— cosh 8 (11) 

with 5 = A + ;A. In equations (9) and (10), A and A are related to 
the resonator geometry and reflectivity parameters by 

sinh A cos A = ^ , (12) 
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cosh A sin A = (2 gi g 2 - 1). (13) 

Furthermore 

72 = «„ (14) 

(cos A + 7 sinh A) /t _ N 

«■-■ = My* ■ — ;HH — rr- (15) 

sin A + cosh A 
The associated eigenvalues K\ {) and K ^ 2> are: 

K? = r exp (-iM)^ [ 1 + , '^ ( ~^ J ' exp (-n«), (16) 

«« = exp (_» f [—flLj]"' exp (-n«). (17) 

The eigenvalue K t which gives the decrease of the reflected field 
after a double pass is the product of the above two eigenvalues given 
by equations (16) and (17). 

K t = K\ X) K\ 2) = (-l)T 1 exp (-j2kd)T Q exp [-(2n + I + 1)5]. (18) 

The eigenvalues K t are exponentially decreasing with I, n. The largest 
eigenvalue is obtained for n = t = 0, corresponding to the fundamental 
TEMo.o mode. The next eigenvalue corresponds to the TEM ll0 mode 
(t = 1, n = 0). Since the eigenvalues are related to the power loss, the 
fundamental mode selectivity will depend primarily on the eigenvalues 
for the TE ,o and TE, l0 modes. 

It is of interest to examine the special case g 2 = 1, which corresponds 
to a resonator with a perfectly reflecting planer mirror M2. This reso- 
nator is completely equivalent to a symmetrical resonator consisting 
of two identical VRM separated by 2d. Both the eigenvalues and the 
fields are the same with the fields beyond d being equal to the reflected 
fields of the nonsymmetrical resonator. 

The modes of the nonsymmetrical resonator are orthogonal at the 
uniform mirror M2, since y 2 = a 2 . Howerer, at the VRM neither the 
incident modes nor the reflected modes are orthogonal. This is shown in 
Appendix A. In addition it is shown that for any particular mode, the 
ratio of the reflected to the incident power at the VRM is precisely 
equal to the absolute value square of the eigen value K t . Physically this 
condition corresponds to conservation of power. 

2.2 Stability Criteria 

For a resonator made to be stable, it is necessary that the exponential 
factors 7 X and y 2 in equations (7) and (8) be finite and have a positive 



2354 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1970 

real part. Both these factors are dependent on cos A and g 2 . The limits 
of the stability regions are thus: (a) g 2 = and (b) cos A = 0. The 
second condition can be expressed in terms of the resonator parameters 
by using equation (13). 

( 2g l g a - 1' 
cash A 



cos A = 1 — I - 



(19) 



(20) 



and the second limit of the stability region is 

1 ± cosh A 

g * = — ^i 

Equation (20) contains the special case of the uniform reflectivity 
resonator (cosh A = 1). For this special case equation (20) reduces to the 
stability criterion derived by G. D. Boyd and H. Kogelnik. 8 

In Fig. 2 illustrative stability diagrams are shown as a function of 
0! and g 2 with exp (2 A) as a parameter. (The choice of this parameter is 
discussed subsequently.) 

A few special cases are considered 

(i) gi = 0. 

This resonator is stable for all values of g 2 , except g 2 = 0. 
(m) g 2 = 0. 

This resonator is unstable independent of the curvature of Ml. 
(Hi) 0i = g 2 = 0. 

This is the very special case of the confocal resonator and is in 

general unstable, except for a URM resonator (fi = 0). 
(iv) 0i = g 2 = 1. 

This is the Fabry-Perot resonator and is stable with a variable re- 
flectivity mirror. 

2.3 The Threshold-Gain Ratio 

To sustain oscillations in a laser resonator a necessary condition is 
that the active medium should have enough gain such that after a 
double transit the field has the same amplitude. This condition can be 
written in terms of the eigenvalues of the resonator modes as 

G \K t \* = 1 (21) 

where G is the power gain per double transit. 

In particular for the TEM 0|0 and TM,, modes, equation (21) can be 
written using equation (18) as: 

G , r^exp(-2A) = 1, (22) 
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Fig. 2— Stability diagrams. 

(? li0 r5exp(-4A) = 1, (23) 

where G ,o and Gi,o is the threshold-power gain required for oscilla- 
tion in the respective modes. A quantity of interest is the threshold- 
gain ratio, t defined by: 



t = j^= exp(2A). 



(24) 



This ratio is a measure of the gain tolerance required for oscillation 
in the dominant TEMo.o mode, and is independent of r . 

The threshold-gain ratio may also be expressed in terms of the loss 
per round trip L n ,o for the TEMn.o mode. Since 



Lo,o = I - Vl exp (-2A). 



(25) 
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The threshold-gain ratio can be written, 



t = 



V 2 

1 o 



1 L .o 



(26) 



Equation (26) is shown in Fig. 3 as a function of L ,o • It is evident 
that the threshold-gain ratio increases with the loss and hence better 
mode discrimination is obtained as the loss increases. 1 Furthermore, the 
power output is related to the power loss per transit. Therefore different 
values of V 2 can be used to shape the spatial distribution of the power 
output. 

A comparison is made (similar to that in Ref . 1) between the threshold- 
gain ratio of a Fabry-Perot resonator with URM and a resonator with 
VRM as a function of the loss per transit. 

Based on the Vainshtein resonator theory 10 for the URM resonator 
the threshold gain ratio 



t = 



1 - L, 



./»... j'-n 



(27) 



where L d is the diffraction loss for the TE ,o mode and i» M>0 is the first 
nonzero root of J„{v) = 0. 

Equation (27) is also plotted in Fig. 3. It is evident from this figure 
that the threshold-gain ratio as a function of the diffraction loss is 
higher for the Fabry-Perot resonator with uniform mirrors than the 
corresponding ratio for the variable reflectivity resonator. This con- 
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Fig. 3— Comparison of uniform and variable reflectivity mirror resonators. 
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elusion was previously reached by Vlasov and Talanov, 1 who also 
showed that the highest threshold-gain ratio is obtained with confocal 
resonators with uniform reflectivity mirrors. However, the high thres- 
hold-gain ratio is only of primary importance for low-gain lasers, 
where the output power obtained by partial transmittivity of the 
mirrors is only a fraction of the power lost by diffraction. For high- 
gain lasers, however, the mode utilization volume is of prime impor- 
tance and the threshold-gain ratio can be kept at a specified level by 
the proper choice of the loss per pass. For such lasers the resonators 
with variable reflectivity mirrors have the advantage that the power 
loss which is necessary for mode discrimination can also be utilized as 
the power output. The mode volume utilization aspect is discussed later. 

III. COMPUTED TBMo.o MODE CHARACTERISTICS 

3.1 Spot Sizes 

The TEMo.n mode is of particular interest since it is the funda- 
mental mode having the highest eigenvalue and hence the lowest loss. 
For this mode the field distributions are gaussian with quadratic phase 
variations. Specifically the field distributions for the TE ,o mode 
from equations (7) and (8) are 

F„ = exp |-~ [exp (-A) cos A + jsinh A sin A]p\\ , (28) 

F rl = exp S -— [exp (A) cos A + j sinh A sin A] P 2 1 > , (29) 

/ Mg 2 / cos A + jsinh A \ 2 \ . . 

F r2 = exp |-— [ s . n A + cogh A J P2 | , (30) 

where F n and Fn are the field distributions of the incident and reflected 
fields at Ml, and F r2 is the reflected field at M2. 

The reflection coefficient can be expressed in terms of A and A is 
by using equation (12) as 

T = r exp (-M/2g 2 sinh A cos A p?). (31) 

The eigenvalue for the TEM .o mode is 

IU = r exp (-j2kd) exp [-(A + fa)] (32) 



with 



^ = J - cos A. (33) 
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The amplitudes of the field distributions at the mirrors are com- 
pletely characterized by the spot sizes defined by that radius when the 
above quantities assume the value of 1/e. Since the exponents in equa- 
tions (28) through (31) are proportional to M, it is convenient to intro- 
duce the Fresnel numbers of the spot sizes. For example, for Fa the 
spot size is defined by 



K 

4#2 



or 



exp (—A) cos Apu = 1 



ir cos A 



(34) 



(35) 



where N iX = p«/Xd is the Fresnel number. The corresponding Fresnel 
numbers of equations (29) through (31), N rl , N r2 , N m are defined 
in an analogous manner. 

The above Fresnel numbers have been computed as a function of the 
threshold-gain ratio t = exp (2A) , with the radii of curvature as param- 
eters. Two types of resonators were considered: (i) a resonator with 
a uniformly reflecting plane mirror and a curved mirror with variable 
reflectivity with radius of curvature as a parameter, and (ii) a resona- 
tor with plane mirror with a variable reflectivity and a uniformly 
reflecting mirror with radius of curvature as a parameter. 

For resonator (i) , Figs. 4 through 6 show the Fresnel numbers of: 
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Fig. 4 — Spot size of incident beam, Nu . 
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Fig. 5 — Spot size of reflected beam, N r 2 . 
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Fig. 6 — Spot size of variable reflectivity mirror, N m . 
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the incident beam at Ml, JV n , the reflected beam at M2, N r2 and the 
Fresnel number of the variable reflectivity mirror N m . Figure 7 shows 
the phase of the eigenvalue ^ , equation (33). Figures 8 through 10 
show the corresponding quantities for the type (u) resonator. The 
phase Vo is the same as in Fig. 7 but with g x and g 2 interchanged. 

A comparison of the characteristics for the two types of resonators 
shows that the most pronounced differences are when either of the 
mirrors have curvatures g x or g 2 = 0.5. Larger spot sizes are obtain- 
able with the type (i) resonator. The Fabry-Pcrot resonator for which 
0i and g 2 = 1.0 is a special case for both types. It also may be noted 
that a large increase in the spot size occurs when one of the mirrors is 
slightly convex, e.g., gr a or g 2 = 1.01. This increase is caused by the 
curvature of the resonator mirror which approaches the unstable 
region, Fig. 2. 

For a finite resonator, the resonator diameter will be limited by the 
minimum obtainable reflectivity at the mirror edges. At the spot-size 
diameter the reflection coefficient of the mirror has the value 1/e. The 
spot-size diameter may be considered as a measure for the diameter of 
the VRM. The Fresnel number of the spot size of the incident beam, 
which for a variety of geometries is the maximum spot size of the 
beam along the resonator is related to the Fresnel number of the 
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Fig. 10— Spot size of variable reflectivity mirror, iV* 
mirror spot size by : 



N r 



= [exp(2A) - 1] = t- 1. 



(36) 



For a resonator with finite dimensions to be a good approximation to 
the infinite resonators, it is necessary that beam power outside the 
mirrors be small. To obtain an estimate of this power, a resonator is 
assumed with a diameter equal to the mirror spot size diameter. 

The ratio p, of the incident power outside the mirror spot-size diam- 
eter to the total incident power is from equations (28) and (1) given 
by 



/ _ exp [-M/2g 2 exp (—A) cos kp\~\pi dp 1 
v = ^£ (37) 

/ exp [— M/2g 2 exp (—A) cos Ap]]pi dpi 

•JO 

After performing the integration and substituting equation (12), 
this ratio can be written as 



V = exp {-2/[exp(2A) - 1]}. 



(38) 
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It readily follows from equation (38) that for a threshold-gain ratio 
t = exp(2A) smaller than 1.43, p is less than one percent. 

A comparison is made between the characteristics of Fabry-Perot 
resonators with one large (such that the diffraction loss is negligible) 
uniformly reflecting mirror and with the other mirror being either of 
uniform or variable reflectivity. For the resonator with VRM, the 
Fresnel number for a given diffraction loss has twice the value than 
that if both mirrors are of the same size. Figure 11 shows the Fresnel 
number of the uniform mirror as a function of threshold-gain ratio. 
The curve is based on the Vainshtein resonator theory. 10 In the same 
figure is also shown the Fresnel number of the spot size of the incident 
field Na • It is evident from this figure that the spot-size diameter at 
the VRM is considerably larger than the diameter of the uniform 
reflectivity mirror for the same values of t. As an example the special 
case of a resonator with a uniform mirror with a Fresnel number of 
two is considered. For this resonator, the field at the mirror has been 
computed by T. Li. 5 Though the field distribution for the TEM ,o is 
not gaussian, for comparison purposes the Fresnel number of the spot 
size (based on the 1/e value for the field) is estimated to be about 
1.6. For the same value of t the Fresnel number of spot size for the 
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Fig. 11 — Comparison of the beam sizes of Fabry-Perot resonators. 
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VRM resonator is 5.2. For lower values of t the difference in the spot 
sizes is even more pronounced. One of the advantages of VRM reso- 
nator is therefore the larger spot sizes and hence the large potential 
for mode volume utilization. 

3.2 Field Distributions in the Resonator 

For efficient mode volume utilization, the field along the resonator 
should be reasonably uniform. The uniformity of the fields is strongly 
dependent on the mirror curvatures. Referring to Fig. 1, let B\ be the 
reflected beam from Ml and B2 the reflected beam from M2. The 
functional dependence of the two beams on the longitudinal z coordi- 
nate has been obtained from the fields at the mirrors, and is given by 
the following equations 



B\ = exp[- 7l (z)]f , 

2 

B2 = exp [-y 2 (z)] | , 



(39) 
(40) 



with 



M\2g 2 \-j exp (A) cos A + j - 



7i(«) = 



(exp (A) + a sin A) : 



+ a 2 cos 2 A — 20 2 - (exp (A) sin A + a) 

z 



(exp (A) + a sin A) 2 + a 2 cos 2 A 



a = 1 + 2gA- - 1 J , 



(41) 
(42) 



and 



M 



2g-. 



~) 2 exp (A) cos A + j[jzr;)\(^ A''' "| A; 



+ cos 2 A + 



d 



72(2) = 



d — z 



[exp (2A)6 + 2 exp (A) sin A(g 2 - 1) - 1] 



exp (A)(& + 



d - z 



+ 



+ 



d -z 

d 
d-z 



- ljsin A 



- 1 cos 2 A 



1] 



(43) 
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with 



6 = (2y 2 - 1). 



(44) 



The real part of equations (41) and (43) has been evaluated as a 
function of z, in terms of the spot-size Fresnel number N x {z) and 
N 2 (z) defined in accordance with equation (34) as N lt o(z) = 2/ 
7i, 2 (z) Ad. For a resonator with g-> = 1.0, Fig. 12 shows the spot-size 
Fresnel number as a function z/d for a number of parameters. 

It is characteristic of resonators with VRM, that minimum-beam 
spot size even for symmetrical resonators does not occur at half the 
mirror spacings in contrast to resonators with uniform reflectivity 
mirrors. In Fig. 12 this characteristic is particularly evident for the 
equivalent confocal resonator g x = 0.5. 

The uniformity of the beams along the resonator increases with the 
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Fig. 12 — Spot size of the beams along the resonator. 
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increasing value of g x up to g x = 1.0 which corresponds to the Fabry- 
Perot resonators. For this resonator the beams are most uniform. For 
the higher value of Qx(Qi = 1.01) the uniformity of the beams within 
the resonator decreases. 

3.3 Minimum Spot Sizes 

The uniformity of the beams within the resonator depends on the 
locations of the minimum spot-size diameters. The beam diameters are 
more uniform if the virtual minimum spot-size diameters occur at 
distances, away from the mirrors, which are large in comparison to 
the separation of the mirrors. The positions of the minimum spot sizes 
are obtained either by determining the maxima or by setting the imag- 
inary parts of equations (41) and (43) equal to zero. Either condition 
gives the same result (i.e., at the minimum spot-size positions the 
beams have constant phase). 

The minimum positions for the two beams {zi/d) mirt and (z2/d) min 
are: 

9 [2 g2 - 1 - exp (A) sin A] , » 

Z?2 [exp (2A) + 2(1 - 2g 2 ) exp (A) sin A + (1 - 2g 2 f] K ' 
and 

(z 2 \ _ , n rexp(A)(2 g2 -l)-sinA] 

ULu " 92 GXP ( } [exp (2A)(1 - 2g 2 ) 2 + 2(1 - 2g 2 ) exp (A) sin A + 1] 

(46) 
The Fresnel numbers of the minimum spot sizes [Ni(z)] min and 
[N 2 {z) ]min are: 

rAr , y, 2g 2 cosAexp(A) ,.~ 

L^WJ-i. - T[exp ( 2A) + 2(1 - 2g 2 ) exp (A) sin A + (1 - 2g 2 ) 2 ] K } 

and 

2 g 2 cos A exp (A) ,.„-, 

[N 2 (z)U» = 7r[ exp(2A)(l-2^) 2 + 2(1 - 2g 2 ) exp (A) sin A + 1] " **> 

The minimum positions and minimum spot sizes have been computed 
for a resonator with a VRM and variable radius of curvature and a 
uniform plane mirror g 2 = 1.0. Since this resonator is equivalent to a 
symmetrical resonator with two VRM, the two minimum positions 
are mirror images with respect to the uniform mirror, and the minimum 
spot sizes are the same. For this resonator 

1 ccvi A 

MM]... = WMU. = - [C08h i _ sin A] • m 
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For specified A equation (49) has a maximum for cosh A sin A = 1, 
which from equation (13) corresponds to the Fabry-Perot resonator 
(g t s g 2 = 1.0). The corresponding Fresnel number N u is: 



N M = 



1 



Tsinh A 



(50) 



Figure 13 shows {z 2 /d) min as a function of g x with t = exp(2A) as 
a parameter. As g x increases, so does {z 2 /d) mm assuming relatively 
large values in the vicinity of g x = 1.0. The large values of (z 2 /d)min 
explain the uniformity of the beams in the Fabry-Perot resonator. 
Figure 14 shows the dependence of the minimum spot on the mirror 
curvature g\ with £ as a parameter. 

3.4 Dependence of the Spot Sizes on the Curvatures of Spherical Mirrors 
The previous calculations show that large spot sizes are obtainable 

in the vicinity of the instability region. How critically the spot sizes 

depend on mirror curvatures is of importance. 
The Fresnel number of the spot size for the incident beam N a 

is directly related to the Fresnel number of the VRM by equation (36) . 




0.6 0.7 0.8 0.9 

CURVATURE OF Ml , g = (l-d/Ri) 



Fig. 13 — Location of the minimum spot size. 
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0.7 0.8 0.9 1.0 

CURVATURE OFMl,g = (l-d/R,) 



Fig. 14 — Fresno] number of minimum spot size as a fund ion of curvature of 
Ml. 

Equations (12) and (13) give the relation between the Fresnel number 
of VRM, the mirror curvature and A. Solving these equations for g-i 
gives : 



r/i = 



ctnh A[(TrN m sinh A) 2 - glf + nN,, 
2wg 2 N m 



(51) 



where N. m is related to /? in equation (1) by N„, = (l/(3\d). 

Equation (51) has been computed as a function of g 2 with Us a 
parameter. Figures 15 through 18 show the computed characteristics 
for N m - 10,20,40,100. 

The critical dependence of the beam spot size N a on the mirror 
curvatures is evident from these figures, particularly as N m in- 
creases. A small change in gi or g 2 results in a large change in 
t and there is a very large change in the beam spot size N a for a 
specified N m . 

The conclusion based on these computations is that though large 
Fresnel numbers for the beam spot sizes are in principle possible, the 
critical tolerance requirements for the mirror curvatures may set a 
practical limit on spot sizes relative to those obtainable with a Fabry- 
Perot resonator with one VRM. 



OPTICAL RESONATORS 



2369 



3.5 Resonator Design 

As an application, a resonator design is considered for a C0 2 laser. 
In view of the realizable high gain per pass, the power loss per pass 
and the related power output should be large. A Fabry-Perot resonator 
with a VRM mirror seems to be most suitable for this application. The 
remaining parameter to be specified is the threshold-gain ratio, t. 
This ratio should be as low as possible in order to obtain large beam 
diameters (see Fig. 11). For the fundamental TEM , o mode operation, 
the limitation on t is based on the accuracy with which the gain can 
be controlled. A value of t of 1,2 is assumed. 

Using the above parameters in equations (28) and (31) (Figs. 4 
and 6) the Fresnel number of the spot size of the VRM is 38.35 and 
that of the incident beam at the VRM is 7.65. For a resonator with 
length d = 100 cm and for a wavelength A = 10 -3 cm (10 micron), 
the radius of the spot size of the VRM, a,,, = 1.95 cm. 

Some of the characteristics of this resonator are shown in Fig. 19. 
Illustrated is the dependence of the power-reflection coefficient T 2 as a 
function of the normalized radius p/a m . Also shown is the normalized- 
incident power density at the VRM, and power loss density at the VRM 
for different values of the reflectivity at the center, T 2 , , as functions of 




1.00 



1.02 1.04 1.06 1.08 1.10 

CURVATURE OF M2, g 2 = (l -d/R 2 ) 



Fig. 15 — Relation between mirror curvature parameters g\ and g s and Fresnel 
number 2V« of spot size at Ml. 
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1.00 



1.02 1.04 1.06 1.08 1.10 

CURVATURE OF M2, g 2 = (l-d/R 2 1 



1.12 



Fig. 16 — Relation between mirror curvature parameters gi and g 2 and Fresnel 
number Nn of spot size at Ml. 




1.02 1.04 1.06 1.08 1.10 

CURVATURE OF M2, g 2 = ( 1 -d/R 2 ) 



Fig. 17— Relation between mirror curvature parameters gi and 17a and Fresnel 
number N tl of spot size at Ml. 
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Fig. 18 — Relation between mirror curvature parameters ffi and g a and Fresnel 
number Nn of spot size at Ml. 

p/a m . The power-loss density is the laser-power output when the 
absorptivity of the VRM is zero. Figure 20 shows the ratio of the power 
loss to the incident power as a function p/a m with Tl as parameters. 

The actual diameter of the VRM can presently be determined only 
by assuming that a finite resonator will behave similarly to a resonator 
with infinite mirrors when the beam power outside a certain diameter is 
small. For the resonator considered, 0.5 percent of the incident beam 
power is contained outside the mirror radius of 0.73 a m and one percent 
outside the radius 0.68 a m which corresponds for the above value of a m to 
1.41 cm and 1.31 cm radii. For a resonator with a VRM of radius a m 
the perturbation of the fields should therefore be very small. 



rv. CONCLUSIONS 

The characteristics of optical resonators with gaussian radial varia- 
tions of the mirror reflectivities have been investigated. These variable 
reflectivity mirror (VRM) resonators seem to be particularly suitable 
for high-gain and high-power laser application such as the 10.6 micron 
C0 2 laser. For the fundamental TEM , o mode generation, these 
resonators have the advantage in comparison to conventional resonators 
that larger beam spot sizes are obtainable (with better mode-volume 



2372 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1970 



t-1.2 

g 1= g 2 = i.o 



POWER REFLECTION 
COEFFICIENT, T 2 




TRANSMITTED BEAM 
_— .^ -DENSITY FOR NO 
' MIRROR ABSORPTION 



0.2 0.4 0.6 0.8 1.0 

NORMALIZED RADIAL COORDINATE, p/a m 



1.2 



Fig. 19— Power distribution at the variable reflectivity mirror. 

utilization) and the power loss necessary for mode discrimination can 
be utilized as the power output. 

The factors limiting the spot are the threshold-gain ratio and the 
mirror-curvature tolerances. 

The Fabry-Perot resonator with a VRM is stable and furthermore 
the field distribution along the resonator is more uniform in diameter 
relative to other resonator geometries. A specific design of such a 
resonator with a gain threshold ratio (Cn.o/Go.o) of 1.2 shows that 
a spot size Fresnel number of 7.65 with a power loss (or power output 
depending on the absorptivity of the mirrors) as high as 40 percent 
of the incident power are obtainable. 

In this investigation it was assumed that the mirrors are infinite. 
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The results presented should be a good approximation to a finite 
resonator when the beam power outside a certain circular region 
has a negligible value (say less than one percent) . 
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APPENDIX A 



Integrals of Laguerre Functions 

In order to determine the orthogonality and power relations for the 
modes in resonators with variable reflectivity mirrors, the following 
integral I' m „ of product of Laguerre functions is evaluated. 

/,< .„ = 2 f exp (-V)U(*P 2 )^(/3pV' +1 dp, 

Jo 

= f X exp (-sOLlMOLimf dt. (52) 

•'0 

For the special case m = n the integral is known. 111 - The integral 
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Fig. 20 — Ratio uf the loss to the incident beam power as a function of p/a m . 
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(52) is evaluated by considering this integral as a Laplace transform 
of two functions f t (t) and f 2 (t) and using the Faltung relation 

r exp (-sO/i(0/ 8 (0 dt = ~ f F 1 ( 2 )F 2 (s - s) dz (53) 

Jo """J Jy—iaa 

where Fi(z) and F 2 («) are the Laplace transforms of f\{t) and f 2 (t), 
y is a constant with Re(s) > Re (7) > 0. 
Let 

U(t) = L!M)= ti-iYQ-D^f (54) 

and 

/ 2 (0 = LtW)t l - (55) 

The Laplace transform of equations (54) and (55) are readily ob- 
tained. Furthermore with the transformation £ = 1/z together with 
equation (53), equation (52) reduces to: 



\7* \ J 2*3? L _iV" +l[(s 



(56) 

In equation (56) the contour of integration encloses the point 
£ = 1/s. Equation (56) is therefore evaluated by determining the 
residue at z = 1/s, which yields 

lL, n = ^7^ (" -)' & {(1 - (*)*[(* - «)r - ir + 'lr-i/.. (57) 
n\ s \ al t/f 

Equation (56) can be expressed in terms of Jacobi polynomials 12 by 
rotating and translating the coordinate system with the result that 



ft (Wjf_0-' /•„ a \m fo _ fftn-mpt.n- 

with 



it.* - {n t +2: (s-a- ens - «— . p:-~<,) (58) 

n! s 



_ s 2 + 2a& - s(a + g) 
' s(s - a - j8) 

and P'„; n ~ m is a Jacobi ploynomial defined by 12 



(59) 



^'^ - 2=fr (1 " * r(1 + * rl 1= [(1 ~ • T)Q+m(1 + x)b+m] - m 
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As an application let F£(p) and F l n *( P ) designate the reflected fields at 
the VRM given by equation (7) and the* indicate the complex con- 
jugate. The integrals which enter in the evaluation of the total reflected 
power due to several modes of the same index I, can be written as 



I 



FUp)FL*(p) P d P = {JM ^ (It. n) r (61) 



where «i is given by equation (11). Using equations (7), (9), (11), (12) 
and (58), it follows that 

fa-*)"' a: „), = t^m+-Mjk exp [- a(« + « + 1 + m 

2 ^••■" n \M cos A * 

I" sin^vl' TsjnM ( _ iA) ]""V:--(,) (62) 
L cos A J L cos A J 

with 

„ = r i + 2 s -Hi^l. (63) 

' L cos A J 

For stable resonators with VRM equation (62) is not equal to zero for 
m ?± n. Hence, the total reflected power is in general not equal to the 
sum of the powers of the individual modes. 

The reflected field F { J } (p) for any particular mode is related to the 
incident field PS ) (p)t by reflection coefficient (1). The evaluation of 
the corresponding integral (61) for the incident fields yields the same 
value for 77. Furthermore setting m = n the following relation if obtained 

ULlz = rj exp [-2Ao(2n + t + 1)] - | K t f '. (64) 

\1 n . n ) i 

The meaning of equation (64) is that the ratio of the reflected to 
incident power for a particular mode is precisely equal to absolute value 
of the eigenvalue squared. 
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